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Abstract 

We present a Monte Carlo algorithm for Hamiltonicity detection in an n-vertex undi- 
rected graph running in O* (1.657") time. To the best of our knowledge, this is the first 
' superpolynomial improvement on the worst case runtime for the problem since the 0*{2^) 

, bound established for TSP almost fifty years ago (Bellman 1962, Held and Karp 1962). 

^jTj' It answers in part the first open problem in Woeginger's 2003 survey on exact algorithms 

^ , for NP-hard problems. 

' For bipartite graphs, we improve the bound to O* (1.414") time. Both the bipartite 

and the general algorithm can be implemented to use space polynomial in n. 

We combine several recently resurrected ideas to get the results. Our main techni- 
cal contribution is a new reduction inspired by the algebraic sieving method for fc-Path 
C/^ ■ (Koutis ICALP 2008, Wilhams IPL 2009). We introduce the Labeled Cycle Cover Sum in 

I which we are set to count weighted arc labeled cycle covers over a finite field of characteris- 

^ • tic two. We reduce Hamiltonicity to Labeled Cycle Cover Sum and apply the determinant 

Q I summation technique for Exact Set Covers (Bjorklund STACS 2010) to evaluate it. 

1 Introduction 

. An undirected graph G = {V,E) on n vertices is said to be Hamiltonian if it has a Hamil- 

I tonian cycle, a vertex order {vo,vi, ■ ■ ■ ,Vn-i) such that fifi+i G E for all i. The indices are 

' enumerated modulo n requiring also that Vn-iVQ is an edge. The problem of detecting if a 

■ graph is Hamiltonian is called the Hamiltonicity problem and was one of the first identified 

, as NP-hard. It is on Karp's original list [16], but is perhaps best known as a special case 

of the Traveling Salesman Problem (TSP). The TSP asks for a tour visiting every vertex 
of an edge weighted graph exactly once that minimizes the total weight. Bellman [3j ^^ and 
^ , independently Held and Karp ^Tj described in the early 1960's a dynamic programming re- 

I currence that solves the general TSP in 0(n^2") time. Their bound also holds for the special 

case of Hamiltonicity and is still the strongest known. Under the widely acknowledged Ex- 
ponential Time Hypothesis, the Hamiltonicity problem has exp{Q{n)) runtime \T^■ There 
is however no known reason to expect the exponential base to be precisely two. Woegin- 
ger in his survey on exact algorithms for NP-Hard problems |22j observes this and asks in 
Open problem 3.1 for a 0*(c") time algorithm for TSP and Hamiltonicity for some c < 2. 
0*(/(n)) suppresses polylogarithmic functions in f{n). We solve the latter problem. 

Theorem 1 There is a Monte Carlo algorithm detecting whether an undirected graph on n 
vertices is Hamiltonian or not running in 0*(1.657"') time, with no false positives and false 
negatives with probability exponentially small in n. 

For graphs having an induced subgraph with many disconnected components, most no- 
tably bipartite graphs, we get an even stronger bound. 
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Theorem 2 There is a Monte Carlo algorithm detecting whether an undirected graph on n 
vertices with a given independent set of size i is Hamiltonian or not running in 0*(2"^*) 
time, with no false positives and false negatives with probability exponentially small in n. 

We also note that our algorithm can be used to solve TSP with integer weights via self- 
reducibility at the cost of a runtime blow-up by roughly a factor of the sum of all edges' 
weights. 

Theorem 3 There is a Monte Carlo algorithm finding the weight of the lightest TSP tour in 
a positive integer edge weighted graph on n vertices in 0*{wl.657"') time, where w is the sum 
of all weights, with error probability exponentially small in n. 

1.1 Previous Work 

Bellman's [3l H] and Held and Karp's algorithm for TSP in an n-vertex complete graph 
G = {V,E) with edge weights i : E ^ M"*" is based on defining Ws^ti^) for s,t £ X Q V as 
the weight of the lightest path from s to t in the induced graph G[X] visiting all vertices in 
X exactly once. This quantity obeys the simple recursion 

o„ { Y\ /™™«GX\{s,t} Ws,u{X \ {t}) + £{ut) : |X| > 2 
"^'^'^^^-[eist] :|X|=2 

Using bottom-up dynamic programming with s fixed, the lightest tour can be evaluated 
by min^gyyi^} uis,t{V) + ^ist) in total 0(n^2") time. An Hamiltonicity instance G can 
naturally be embedded in a TSP instance on the same number of vertices. Simply let the 
weight function i take the value for vertex pairs corresponding to an edge in G, and 1 
otherwise. 

Another algorithm amenable to Hamiltonicity with (almost) the same running time 
is the inclusion-exclusion counting over n-long closed walks in the induced subgraphs. The 
algorithm has been (re) discovered at least three times \18\ 117 ^ [2], Let s be any vertex in the 
graph, then the number of Hamiltonian cycles is given by 

E (_i)in(xuM)i(A[xuM]").,. 

XCV\{s} 

Here, A[Y] denotes the adjacency matrix of the induced graph G[Y], and (A[y]")s^s the entry 
at row and column s of the matrix A[y]". The idea behind the algorithm is that crossing 
walks will be canceled since they are counted equally many times with the sign factor +1 as 
with —1. 

In restricted graph classes the general 0(2") bound has been sharpened. Broersma et 
al. [8] proved that Hamiltonicity in claw-free graphs has an O* (1.682") time algorithm. 
Iwama and Nakashima |15) improving slightly on Eppstein [I2|, showed that TSP in cubic 
graphs admits an O* (1.251") time algorithm. In graphs of maximum degree four, Gebauer |13] 
described how to count the Hamiltonian cycles in 0*(1.715") time. For larger degrees, only 
minuscule improvements are known. Bjorklund et al. [7] observed that both the Bellman- 
Held-Karp dynamic programming and the inclusion-exclusion algorithm need to look only at 
X for which G[X] is connected and include s. These are at most (2 — e)" with e depending 
inversely exponentially on the maximum degree. 
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Another line of research addresses a natural parameterized version of the problem called 
the A;-Path problem: how much time is required to find a simple (noncrossing) path on k 
vertices in an n vertex graph. Alon et al. [1] showed the first c'^n*^^^^ time algorithm for 
some constant c for the problem. The constant c has since been improved several times 
culminating in the work of Koutis [19] introducing an interesting algebraic sieving technique. 
His algorithm was subsequently refined by Williams [21] to yield a 2*^n'^(^) time algorithm 
for the /c-Path problem. In the extreme k = n, their algorithm's runtime coincides with the 
previously best 0*(2") time bound for Hamiltonicity. 

1.2 Our Approach 

The inclusion-exclusion algorithm in the previous section has several desirable properties: It 
uses space polynomial in the input size, it works also for directed graphs, it is deterministic, 
and it is capable of counting the solutions. Our algorithm also manages with only polynomial 
space after some extra care, but does not have the other three properties. Indeed, we crucially 
depend on the graph being undirected and that polynomial identity testing has an efficient 
randomized algorithm (whereas no deterministic is known). Moreover, our strategy does not 
even seem to be able to approximate the number of solutions. 

The inclusion-exclusion algorithm can be thought of as first counting too much (all closed 
n- walks through s) and then canceling out every false contribution (crossing n- walks). We 
will take a similar approach, but count instead weighted cycle covers in directed graphs over 
fields of characteristic two. A cycle cover in a directed n- vertex graph is a set of n arcs such 
that every vertex is the origin of one arc, and the end of another. The arcs together describe 
disjoint cycles covering all vertices of the graph. In particular, the cycle covers contain the 
Hamiltonian cycles. 

Our algorithm is much inspired by the recent work of Koutis [19] and Williams [21] for 
A;-Path. Although we don't apply their work directly, we use several of their ideas. In 
particular we evaluate multivariate polynomials over fields of characteristic two to sieve for 
the Hamiltonian cycles, just as they do for A;-paths. 

The main new technical ingredient is the introduction of determinants to count weighted 
cycle covers. This is an extension of the idea of using determinants to count perfect matchings 
recently employed by Bjorklund [5J for Exact Set Cover. Unlike Koutis-Williams we are 
unable to construct small arithmetic circuits. Instead we depend on the efficient algorithms 
for computing a matrix determinant numerically. 

1.3 Organization 

The rest of the paper is organized as follows. In Section[2]we introduce the technical machinery 
needed. In particular we define our weighted cycle cover problem, hint at how it relates to 
the Hamiltonicity problem, and presents a way to compute it. In Section [3] we describe 
how we can reduce a Hamiltonicity instance to the weighted cycle cover problem and prove 
our main claims Theorem [1] and [2] given exponential space. In Section [J] we argue how the 
algorithm can be modified to use only polynomial space. Finally, in Section [5] we note how 
the technique can be extended in a known way to solve TSP to prove Theorem [3l 
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2 Preliminaries 



In a directed graph D = {V, A) a cycle cover is a subset C ^ A such that for every vertex 
V £ V there is exactly one arc a^i S C starting in v, and exactly one arc 0^2 £ C ending in 
V. The graphs in this paper have no loops, i.e. arcs connecting a vertex to itself, and thus 
we also have that ayi 7^ 0^,2- We denote by cc{D) the family of all cycle covers of D, and by 
hc{D) C cc{D) the set of Hamiltonian cycle covers. A Hamiltonian cycle cover consists of one 
big cycle passing through all vertices. The remaining cycle covers (which have more than one 
cycle), cc{D) \ hc{D), are called non-Hamiltonian cycle covers. Although an element of hc(D) 
is formally a subset of arcs, we will sometimes write it as a vertex order (vo,vi, ■ ■ ■ ,Vn-i) 
implicitly referring to the arcs WjWj+i as the actual Hamiltonian cycle. For undirected graphs 
G, hc{G) includes the Hamiltonian cycles with orientation, i.e. traversed in both directions. 
Hence we will for a Hamiltonian cycle H S hc{G) for an undirected graph G, talk about arcs 
uv £ H inferring that the cycle is oriented from n to f along the edge uv. 

We write g : A ^ B for a surjective function g from the domain A to the codomain 
B. For a function g : A ^ B we associate the function g^^ : B ^ 2^ as its preimage, 
g'~^{b) = {a £ A: g{a) = b}. For a matrix A, we denote by Aij the element at row i 
and column j. For a polynomial p{r) in an indeterminate r we write [r']p(r) to address the 
coefficient of the monomial r' in p{r). 

We will reduce Hamiltonicity to a variant of cycle cover counting defined next. We 
introduce the Labeled Cycle Cover Sum. The name stems from the fact that every arc 
in the cycle cover is labeled by a nonempty subset of a set of labels. 

Definition The Labeled Cycle Cover Sum for a directed graph D = {V,A), a label set 
L, and a function / : A x 2^ \ {0} — )• on some codomain ring R is 

A{D,LJ)= E Uf(^^9-\a)). (1) 

C(^cc{D) g:L^C aeC 

Note in particular that the inner sum is over all surjective functions g, meaning that the 
label g~^{a) is a nonempty subset of L for all arcs a € C. In words the computation is over 
all arc labeled cycle covers of the graph such that all arc labels are nonempty, are pairwise 
disjoint, and together exhaust all of the labels L. 

2.1 Cycle Cover Cancelation in Characteristic Two 

In this section we show that restrictions on the graph, the computation ring, and the function 
/ can be imposed so that the resulting summation in a Labeled Cycle Cover Sum instance 
is over the Hamiltonian cycle covers only. The contributions of the non-Hamiltonian cycle 
covers are canceled out. First, we say a directed graph is bidirected if it for every arc uv 
has an arc in the opposite direction, vu. Second, we let the ring R have characteristic two. 
Third, for an arbitrarily chosen special vertex s, f : A x 2^ \ {0} — )• i? is an s-oriented mirror 
function if f{uv,Z) = f{vu,Z) for all Z and all u ^ s,v ^ s. This definition asymmetry 
around the vertex s is a first precaution to avoid that everything cancels. We still want the 
contributions of the Hamiltonian cycle covers to leave a nonzero result. We will address this 
aspect further in Section 12.21 

The following lemma captures how the non-Hamiltonian cycle covers vanish, which also 
will imply the nonexistence of false positives in our resulting algorithms. 
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Figure 1: Two pairs of labeled non-Hamiltonian cycle covers mapped by M in the proof of 
Lemma m (a) and (b) are dual, the arcs along the cycle C in (a) are reversed in (b) but keep 
the same labeling. The fixed-point free duality works also when the cycle consists of only 
two vertices as shown in the pair (c) and (d). Note also that this would not be the case if the 
labels a and b were allowed to be the empty set, since if so they would coincide. 

Lemma 4 Given a bidirected graph D = (V,A), a finite set L, and special vertex s £V, let 
f be an s-oriented mirror function with a codomain ring of characteristic two. Then 

A{D,L,f)= Yi E n/(«'5"'(«))- 

Hehc{D) giL^H aeH 

Proof Confer the definition of Labeled Cycle Cover Sum ([T]). A labeled cycle cover is a 
tuple (C, g) with C G cc{D) and g : L ^ C. We will argue that the labeled non-Hamiltonian 
cycle covers can be partitioned into dual pairs such that both cycle covers in every pair 
contribute the same term to the sum. Since we are working in a ring of characteristic two, all 
these terms cancel. To this end we define a mapping M from the labeled non-Hamiltonian 
cycle covers onto themselves. 

Consider a labeled non-Hamiltonian cycle cover {C,g). We define M{C,g) = {C',g') as 
follows. Let C be the first cycle of C not passing through s. Note that there must exist one 
since the cycle cover consists of at least two cycles and all cycles are vertex disjoint. Here 
first refers to any fixed order of the cycles. Let C = C except for the cycle C which is 
reversed in C", i.e. every arc uf G C is replaced by the arc in the opposite direction vu in 
C . Note that this arc exists since the graph D is assumed to be bidirected. In the special 
case when C consists of only two arcs, C is identical to C . The function g''^^ is identical to 
g~^ on C \ C, and is defined by g'^^{uv) = g~^{vu) for all arcs uv G C. In other words, the 
reversed arcs preserve their original labeling. Note in particular that {C,g) ^ M{C,g) and 
{C,g) = M{M{C,g)). Hence the mapping M uniquely pairs up the labeled non-Hamiltonian 
cycle covers (cf . Fig. [1]) . 

Since / is an s-oriented mirror function and has f{uv,Z) = f{vu,Z) for all arcs uv not 
incident to s and all Z G 2^ \ {0}, {C,g) and M{C,g) contribute the same product term to 
the sum in ([1]) and hence cancel. 

2.2 Detecting the Hamiltonian Cycles 

In the previous section we argued that the non-Hamiltonian cycle covers' contributions to the 
sum in the LABELED Cycle Cover Sum cancel if certain requirements are met. For this to 
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be useful we also need that the Hamiltonian cycle covers don't vanish. To this end, it will be 
instructive to think of the elements of / as nonconstant multivariate polynomials in variables 
associated with the argument arc and label set. In particular for elements adjacent to the 
special vertex s, f{su,X) and f{us,X) will not share variables for any su,us G A. This will 
ensure that the Hamiltonian cycles oriented in opposite directions will contribute different 
terms to the sum. 

In Section [3?T] and [32] we wih define / such that the associated Labeled Cycle Cover 
Sum seen as a polynomial in the underlying variables will have at least one unique monomial 
per (oriented) Hamiltonian cycle. Moreover, there will be no monomials resulting from non- 
Hamiltonian cycle covers, as a consequence of Lemma [H 

To detect if the polynomial resulting from the Labeled Cycle Cover Sum is identically 
zero (=no Hamiltonian cycles) or not (=at least one Hamiltonian cycle), we will employ the 
old fingerprint idea often attributed to Freivalds (see [2D] for a historical account). We will 
evaluate the polynomial in a randomly chosen point. If the fingerprint result is nonzero we 
know for sure the polynomial couldn't possibly be the zero polynomial. If the result is zero 
we guess that so is the polynomial. The Schwartz-Zippel Lemma (see e.g [201 P- 165]) ensures 
that with great probability we will succeed: 

Lemma 5 (Schwartz-Zippel) Let P{xi,X2, ■■■,Xn) be a nonzero n-variate polynomial of 
total degree d over a field F. Pick ri,r2, ...,r„ G F uniformly at random, then 

Pr(P(ri,r2,...,r„) =0) < ^ 

Note that the algorithms' actual computations in this paper will be over a finite field 
GF{2^) for some k after replacing the variables for values. In the analysis in Section [3. II and 
13.21 though, the function / will be treated as a multivariate polynomial with coefficients from 
GF{2). 

2.3 Determinants and Inclusion— Exclusion 

Bjorklund in [5] presented a computation technique which suitably tuned can be used to 
solve Labeled Cycle Cover Sum relatively quickly. It relies on the well-known fact that 
the determinant of an n x n-matrix A over a ring of characteristic two coincides with the 
permanent. 

n 

det(A) = per(A) = ^ H^^-W (2) 

(7:[n]->[n] i=l 

The summation is over all permutations o" of n elements. 

Permanents have a natural interpretation as the sum of weighted cycle covers in a directed 
graph. Formally let D = {V, A) be a directed graph with weights w : A ^ R, and define a 
\V\ X \ V\ matrix with rows and columns representing the vertices V 




w{ij) :ijeA 
: otherwise 



then 

per(A)= n^(«)- (3) 

CGcc(D) aeC 
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We will see that Labeled Cycle Cover Sum can be evaluated through a sum of an 
exponential number of determinants. To this end we define for every Z C. L the matrices 

''J \ : otherwise. ^ ^ 

We introduce a polynomial in an indeterminate r, with r aimed at controlling the total 
rank of the subsets used as labels in our labeled cycle covers. 

p(/,r)= J]det I J]rl^lM/(Z) I (5) 

YCL \ZCY I 

This polynomial can be thought of as an inclusion-exclusion formula in disguise, which actu- 
ally computes an associated Labeled Cycle Cover Sum in characteristic two. 

Lemma 6 For a directed graph D, a set L of labels, and any f : A x 2^ \ {0} — )• GF{2^) 

[r\%if,r)=AiD,L,f). 

Proof Rewriting the expression for p(f, r) ^ using the equivalence of the determinant and 
the permanent in rings of characteristic two ([2]), the cycle cover interpretation of the perma- 
nent ([3]), and the matrices Mj we get 

YCL Cecc{D) q:C^2Y\{9} aeC 

Changing the order of summation, we have 

C€cc{D) q:C^2l^\{9} \JaeC "GC 
YCL 

For functions q : C ^ 2^ \ {0} such that UaGC ^('^) ^-^^ '^^ose union over the elements 
doesn't cover all of L, we note that the innermost summation is run an even number of times 
with the same term (there are 2^^'^^^ec '^^"'^^ equal terms). Again, since the ring characteristic 
is two, these cancel. We are left with 

pif^^)= E E 

and in particular, the coefficient of rl^l 

[r^%{f,r)= E n 

C<^cc(D) g:C^2^\{0} a&C 
Va^6:g(a)ng(fe)=0 

since lJaGc'?('^) ~ ^ SaeC ~ 1^1 implies Va / 6 : q{a) H q{h) = 0. Inverting the 

function q we arrive at the Labeled Cycle Cover Sum definition 
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The above lemma is the base identity enabhng a relatively efficient algorithm for comput- 
ing Labeled Cycle Cover Sum. The runtime is exponential in the number of labels, but 
polynomial in the size of the input graph. 

Lemma 7 The Labeled Cycle Cover Sum A{D, L, f) for a function f with codomain 
GF{2^) on a directed graph D on n vertices, and with 2^ > \L\n, can be computed in 
0((|-Lpn + |L|n"^+")2'^l + |Lpn^) arithmetic operations over GF{2^), where ui is the square 
matrix multiplication exponent. 

Proof We evaluate the Labeled Cycle Cover Sum via the identity in Lemma El We ob- 
serve that p{f, r) as a polynomial in r has maximum degree \L\n. To recover one of its coeffi- 
cients (the one for r'^'), we need to evaluate the polynomial for \L\n choices of r and use inter- 
polation to solve for the sought coefficient. We can for instance use a generator g of the multi- 
plicative group in GF{2^) and evaluate the polynomial in the points r = g^,g^,g^, . . . . , g^^^^~^ . 
The requirement 2^^ > \L\n ensures the points are distinct, and hence that the interpolation 
is possible. For every fixed r, our algorithm begins by tabulating T{Y) = Ezcy ^'^'^/(^) 
for all y C L through Yates' fast zeta transform |23j (see also for a recent treatment) 
in 0(|L|2l^l) field operations. Next we evaluate p{f,r) = Eycl det(r(y)) in 0(n"2l^l) op- 
erations using the determinant algorithm by Bunch and Hopcroft [9], with uj = 2.376 the 
Coppersmith- Winograd square matrix multiplication exponent [10], and store the value ob- 
tained. Summing up the number of field operations required over all \L\n values of r, the first 
part of the runtime bound follows. Once all values are computed, we employ the 0(|Lpn^) 
time Lagrange interpolation. 

3 The Reduction 

We wiU reduce Hamiltonicity to Labeled Cycle Cover Sum. The overaU idea is to 
partition the vertices of the input undirected graph G into two equal halves. We construct a 
smaller bidirected graph D on one of the halves, and use the other half as labels in a Labeled 
Cycle Cover Sum on D. An especially simple reduction is possible for bipartite graphs, 
which we describe next, even though the result will also follow directly from Theorem [2j 

3.1 Warm-up: Hamiltonicity in Bipartite Graphs 

We are given an undirected bipartite graph G = {Vi,V2,E) on n vertices. We describe an 
0*(2"'/^) time algorithm detecting if G is Hamiltonian. We know a Hamiltonian cycle if it 
exists will alternate vertices from Vi and V2 along the cycle. Thus we can safely assume = 
IV2I = n/2 since otherwise the graph is not Hamiltonian. We will follow the setup outlined 
in the previous section and imagine a symbolic Labeled Cycle Cover Sum describing a 
multivariate polynomial over a ring of characteristic two. We fix a special vertex s £ Vi 
and introduce variables Xuv and x^u for every edge uv G E. We equate 
when u = s or V = s. For every pair of different vertices u,v S Vi, we define the set 
N{u,v) = {w: w G V2,uw £ E,wv £ E}. We set D = {Vi,F) with F including arcs in 
both direction for every pair of different vertices u,v £ Vi such that N{u, v) 7^ 0. For an arc 
uv G F, and a vertex w G N{u,v), we set f{uv,{uj}) = XuwXwv In all other points / is set 
to zero. 

Lemma 8 With G,D,V2, and f defined as above, 
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/ A{D, V2, f) = J2Hehc(G) UuveH ^uv 

II A{D, V2, f) is the zero polynomial if and only if hc{G) = 0. 

Proof I. Since D is bidirected and / is easily seen to be an s-oriented mirror function, we 
have from Lemma H] that 

A(A^2,/)= E \{f{^^9-\a)) 

Hehc{D) g:V2^HaeH 

Since f{a,X) is nonzero only when X is a single vertex in V2, we can rewrite the identity as 

A{D,V2,f)= Yl E n 

H(^hc{D) q:H->-V2 a€H 

Here the summation is over all functions q which are one-to-one. Replacing / by its definition 
we get 

HD,V2,f)= X] X] n ^wq(wz)Xq(wz)z 

H£hc{D) q.H^V2 wzeH 
Vuv(^H: 

Since q is one-to-one, every vertex in V2 is mapped to by precisely one arc G F on a 
Hamiltonian cycle in D, and we have 

A{D,v2,f)= n ^™ 

H&hc{G) uveH 

as claimed (remembering that hc{G) contains all oriented Hamiltonian cycles in an undirected 
graph G). 

H. Clearly from I, A(i^,V2,/) is zero if G is non-Hamiltonian. In the other direction, 
we argue that every undirected Hamiltonian cycle will contribute two different monomials 
each to the sum. This is again because of the special vertex s. Every monomial term in 
the summation corresponds to an oriented Hamiltonian cycle in G and the variables of the 
polynomial are uniquely associated with one edge of the graph. Two monomials resulting 
from two different undirected Hamiltonian cycles will have some variable the other doesn't 
have. Thus the only chance of two monomials being identical would be the pair of monomials 
resulting from the same undirected Hamiltonian cycle in opposite orientations. Since the 
variables tied to the oppositely directed arcs incident to s are different, these are also unique 
monomials in the sum. 

3.1.1 Algorithm and Analysis 

The algorithm repeats the following process, called a run, a number of times linear in n: 

The setup for Lemma [8] shows how to transform the input graph G = {Vi,V2,E) with 
l^il = 1^2! = "'/2 into a symbolic Labeled Cycle Cover Sum A(Z), L, /) on a graph D on 
n/2 vertices and n/2 labels L. We set k large enough, say 2^ > cn for some c > 1. Next we 
evaluate A(D, L, /) in a randomly chosen assignment point p to the variables over the field 
GF{2^) with the algorithm from Lemmad Lemma [5] tells us it will with probability at least 
1 — 1/c result in a nonzero answer if and only if A(Z), L, /) was a nonzero polynomial (and 
G Hamiltonian). If any run results in a nonzero answer, we output that G is Hamiltonian, 
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otherwise not. Since the algorithm uses a hnear number of runs in n the probabihty of false 
negatives is brought down to exp{—^l{n)). 

The time to compute A{D,L,f) in p is dominated by the runtime in Lemma [7] and is in 
0*(2"'/^). Summing over all runs, the total time bound only grows by a factor linear in n. 

3.2 The General Case 

In a general Hamiltonian undirected graph G = {V,E), unlike the bipartite case, we don't 
know a priori which subset of the vertices will be traversed every other vertex along a Hamil- 
tonian cycle in G. Hence it is difficult to partition the vertices in two equal parts of which 
one could serve as labels as in the previous section. Still, a uniformly randomly chosen par- 
tition V = ViL)V2 with I Vi I = I V2 1 has with large enough probability a property that we can 
exploit: The number of transitions along a fixed Hamiltonian cycle from a vertex in one part 
to a vertex in the other part is n/2 in expectation. Note that in the bipartite case it was n. 
The vertices in Vi in the bipartite case were handled at a polynomial time cost whereas the 
ones in V2 came at a price of a factor 2 each in the runtime. In the same vein, the vertices 
in Vi followed by a vertex in V2 along a fixed Hamiltonian cycle in the general case will be 
computationally cheap. To see how, we need to distinguish arcs along the fixed Hamiltonian 
cycle according to the partition. 

For a Hamiltonian cycle H, we call arcs connecting adjacent vertex pairs along 
H unlabeled by V2 if both Vi and Uj+i belong to Vi. The remaining arcs are referred to as 
labeled by V2. We partition the arcs of H in C{H) as the set of labeled, and U{H) as the 
set of unlabeled arcs by V2. We will use that there aren't too many arcs unlabeled by V2. 
Define hcy^ (G) as the subset of hc{G) of Hamiltonian cycles H which have precisely m arcs 
unlabeled by V2 along H. 

We introduce variables x^v and x^u for every edge uv £ E such that u G V2 or v £ V2 (or 
both). We identify Xuv with x^u except when u = s ov v = s. 

We consider a complete bidirected graph D = (Vi, F) and use V2 as some of the labels. In 
addition to V2 we will add a set Lm of size m of extra labels aimed at handling arcs unlabeled 
by V2. For each edge uv in G'[yi] and every element d € we also introduce new variables 
Xuv,d and Xvu,d- Again we let Xyt,^^ coincide with x^^^^ except when u = s or v = s. 

For two vertices u,v G V, and a nonempty subset X <^V we define Vu,vi^) as the family 
of all simple paths in G from u to v passing through exactly the vertices in X (in addition to 
u and v). For uv £ F and C X C V2, we set 

f{uv,X) = ^ a;^^ 
For every arc uv £ F such that uv is an edge in and every d G L^, we set 

f{uV,{d}) = Xuv,d 

In all other points / is set to zero. 

Lemma 9 With G,D,V2,U,C,m,Lm and f defined as above, 

I A{D, V2 U Lm, f) = Y^Hehc-^^{G) (^a:U(H)^Lm IluveU{H) ^uv,a{uv)^ (lluvGCiH) 

with a one-to-one. 
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// A{D, V2 U Lm, f) is the zero polynomial if and only if hcy^ (G) = 0. 

Proof I. Since D is bidirected and / is an s-oriented mirror function, we have from Lemma H] 
that 

Hehc{D) g:V2UL,n^H aGH 

A Hamiltonian cycle H G hc{G) is naturahy associated with a Hamiltonian cycle in D by 
simply omitting the vertices along H which belong to ¥2- We want to go the other way, to 
expand the Hamiltonian cycles in D into Hamiltonian cycles of G. To do this we observe 
that the arcs of a Hamiltonian cycle in D in the sum above either are labeled by an element 
of Lm or a nonempty subset of V2, since these are the only subsets of the labels for which 
/ is nonzero. We extend the definition of labeled and unlabeled arcs (which were defined 
previously for Hamiltonian cycles in G only). For a Hamiltonian cycle H £ hc{D) labeled by 
the function g : V2U Lm H we say an arc uv £ H is labeled by V2 if g^^{uv) C V2, and 
unlabeled by V2 if g^^{uv) G Lm- 

Since every arc unlabeled by V2 along a Hamiltonian cycle consumes exactly one of the 
m labels in Lm, and all labels are used, we have that only the Hamiltonian cycles in D with 
exactly m arcs unlabeled by V2 leave a nonzero contribution. We expand the summation in 
all possible labeled and unlabeled arcs of the Hamiltonian cycles in D. We note that a cycle 
H leaves a nonzero result only if the m arcs unlabeled by V2 along the cycle are also edges in 
G. 

A{D,V2ULmJ)= Y E ^Hu{Lm)AH,{V2) 

H&hc{D) HuUHl=H 



HuDHl-- 
\Hu\=m 



with 



and 



,cr:Hu-^Lm a&Hu 



\g.V2^HL a£HL 

Here the summation is over all functions a which are one-to-one. Replacing / by its definition 
we further expand the inner expressions to 

AHui^m) = j Yl n ^uv,a(uv) 

and 

^hav2)=1 En E n 

Every vertex in V2 is mapped to by precisely one arc in F on a Hamiltonian cycle H m. D. 
Rewriting the expression as a sum of Hamiltonian cycles in G we have 

A{D,V2yjLmJ)= E E n 

H(^hc^^{G) \cT:U{H)^LmUveU{H) / \uveC{H) 
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as claimed. 

II. If the graph G has no Hamiltonian cycles, the sum is clearly zero. For the other 
direction, we see that each Hamiltonian cycle contributes a set of ml different monomials 
per orientation of the cycle (one for each permutation a), in which there are one variable 
per edge along the cycle. Monomials resulting from different Hamiltonian cycles thus are 
different since two different Hamiltonian cycles each has an edge the other has not. Every 
pair of oppositely oriented Hamiltonian cycles along the same undirected Hamiltonian cycle 
also traverse s through oppositely directed arcs. Since the variables tied to the oppositely 
directed arcs incident to s are different, these are also unique monomials in the sum. 

3.2.1 Algorithm 

We are given an n vertex undirected graph G = (Vi E) as input, assuming n even for simplicity. 
The algorithm repeats the following process, called a run, several times r to be specified later 
in Section [3.2.21 

In each run, a partition ViU V2 = 1^ is picked uniformly at random, with |Vi| = = n/2. 
Next the algorithm loops over m, the number of edges unlabeled by V2 along a Hamiltonian 
cycle, from through rrimax, where rumax is specified in Section r3.2.2l For each value of m, 
the setup for Lemma [9] describes how to transform the input graph G given m,Vi,V2 into a 
symbolic Labeled Cycle Cover Sum A(D, V2 U Lm, f) on a bidirected graph D on n/2 
vertices and n/2 + m labels V2 U Lm- We operate over a field GF{2^) with k again set large 
enough, say 2^ > cn for some c > 1. 

Next in each run, a point p assigning values from the field GF{2^) to the variables in / is 
chosen uniformly at random. The function / in p is tabulated for all subsets of V2 (in other 
points the function is easy to compute). This can be achieved by running a variant of the 
Bellman-Held-Karp recursion. Formally, let f : {V xV) x. 2^^ — t- GF(2'^) be defined for u ^ v 
and C X C ^2 by 

f{uv,X) = ^ X^z 

Peru,v{x) wzeP 

then f{uv, X) = f{uv, X) for ut; G F and C X C V2. 
For |X| > 1 the recursion 

f{wv,X\{w}) 

w£X,uw£E 

can be used to tabulate f{.,X) for all C X C V2. 

Finally, in each run, we evaluate A(D, V2 U Lm, f) in p through Lemma[7l If the result in 
any of the runs is nonzero we conclude and output that G is Hamiltonian, otherwise not. 

3.2.2 Analysis 

Since the tabulation for / is in 0*{2^'^^) time, the runtime of each run is dominated by 
the runtime in Lemma [71 The worst case occurs for m = rUmax in which case we get a 
0*(2°-^"+™""''^) time bound. The total runtime is in 0*(r2°-^'^+™'"''^), and the probability of 
false negatives is at most PriYZZo" \hcv^{G)\ = 0)^ 

Lemma [S] tells us that every run will with probability at least 1 — 1/c result in a nonzero 
answer if and only if A(Z), V2 U Lm, f) was a nonzero polynomial (and G has a Hamiltonian 
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cycle with m edges unlabeled by ¥2). A straightforward application of Markov's inequality 
bounding the probability that a fixed Hamiltonian cycle gets more than n/4 edges unlabeled 
by V2 shows that Pr(X;^=o l^c^a^'^)! = 0) ^ n/(4 + n). From this we can deduce that setting 
rrimax = n/A and r = il(n^), wiU give a total runtime bound of 0*(2t'^) 0*(1.682") and 
exponentially small probability of failure in n. 

A slightly better runtime bound is obtained by trading the probability of success in a 
single run for more runs (suggested to the author by Ryan Williams and Petteri Kaski, 
independently of each other). To see just how much better we need a stronger bound on the 
probability. 

Lemma 10 Let G = {V, E) he a Hamiltonian undirected graph, and Vi |J V2 = ^ with |Vi| = 
1^2! = n/2. Then 

Proof We will obtain the first bound by counting the probability that one fixed Hamiltonian 
cycle H = (i;o,fi, • • • , Vn-i) has exactly m arcs unlabeled by V2) and moreover has vq G Vi 
and Vn-i G V2. For such a H there are exactly n/2 — m indices i such that Xi G Vi and 
Xj+i G V2, and just as many indices i where Xi G V2 and Xi^i G Vi. The ordered list of these 
transition indices 11,12,- " > ^n-2m uniquely describes the partition, and in the other direction 
any such list with < ii, Vj : ij < ij+i, and in-2m = n — 1 corresponds to a unique partition. 
Set iQ = —1 and define the positive integers dj = ij — ij-i for all < j < n — 2m. Note that 
di,d3, - ■ ■ , dn-2m-i describes a partition of the vertices in Vi in n/2 — m groups. Analogously, 
d2,d4, - ■ ■ , dn-2m describes a partition of the vertices in V2 in n/2 — m groups. The number 
of ways to write a positive integer p as a sum of k positive integers is (^Zj) • Multiplying the 
number of ways to partition the vertices in Vi with the number of ways for V2, and dividing 
with the total number of balanced partitions (^"2) ' ^'■^^ result follows. 

The second bound is derived by replacing the binomial coefficients with their factorial 
definition and using Stirling's approximation for n! G 0((n/e)"-v/27rn). 

To bring the probability of false negatives down to exp{—Q{n)), we need the number of 
runs r = n^'^^^Pr~^{\hcy^°-'^ {G)\ > 0). Solving for the local minimum of the total runtime, 
using Lemma [To] to bound the probability, we get rrimax = 0.205 and r = 7^0(1)20. 024n ^.^j^g^ 
Altogether, a runtime bound of O* (1.657"). 

For Theorem [21 we set V2 equal to the independent set of size i given as input instead of 
the random partition in the first step of the algorithm above. We note that a Hamiltonian 
cycle must have at least 2i arcs labeled by V2. This is because every vertex of the independent 
set must be incident to two arcs along the cycle, and no arc is connected to more than one 
of them since they are disconnected by definition. Hence, if i > n/2 the graph is surely 
non-Hamiltonian, and if z < n/2 we only need the loop over m in the above algorithm to 
count to rumax = n — 2i, and let r be linear in n. This gives us the 0*(2"~*) time bound. 

4 Polynomial Space 

The algorithm in Lemma[7| invokes Yates' fast zeta transform which uses almost as much space 
as time. The tabulation of the function / also uses space exponential in n. Here we describe 
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an alternative way of solving the problem with an algorithm using only polynomial space. 
This also enables the summation task in the Labeled Cycle Cover Sum to be divided on 
several processors in parallel. We adopt the notation from Section 13.21 There we described 
/ as a function of simple paths in G[V2]- The idea here is to embed the inclusion-exclusion 
counting over walks to sieve for the simple paths. We let yVu,v{X,^) be the set of walks of 
length I in G starting in u and ending in v but in-between visiting only vertices from X. For 
a walk W £ yVu,v{X,l) we define its support on X, denoted SxiW), as the set of vertices in 
X traversed by the walk. With W = {u,wi,W2,---, wi-i,v) we have Sx{W) = IJj=i ^j- 

Note that Vu,v{X) C >Vu_t,(X, |X| + 1). We will derive an analogue g of the function 
/. We replace simple paths for walks. The point being that the inner summation in p(f, r) 
in ([5]) will take the form of counting over walks, which we know how to do fast using only 
polynomial space. 

For an arc uv G F, and a subset C X C V2 we set 

n 

g{uv,X,r)= ^ ^ JJ x^, 

k=\X\ W&Wu,viX,k+l)wz€W 
Sx(W)=X 

For every arc uv £ F such that uv is an edge in G[Vi], and every d G Lm, we set 

g{uv,{d},r) = rXuv,d 

In all other points g is set to zero. 

Remembering the definition of M from we define an analogue of r) from ([5]). 

g(g,r) = J^detl J]M,(„„,)(Z)| (6) 

Lemma 11 For an undirected graph G = {V,E), a vertex partition ViL)V2 = V, and with f 
defined in Section l^.i^l and g defined as above, it holds that 

[r^]q{g,r) = [r^]p{f,r) 

Proof Rewriting the expression for q{g, r) ([6]) using the equivalence of the determinant and 
the permanent in rings of characteristic two ([2]), the cycle cover interpretation of the perma- 
nent ([3]), and the matrices Mg( ^ Q, we get 

(l{g.r)=Y^ ^ ^ J\g{a,h{a),r) 

YCL C€cc(G) h:C^2^\{0} a£C 

Changing the order of summation, we have 

CecciG) h:C^2L\{lli} U^g^ hia)CYCL a£C 

For functions h : C ^ 2^ \ {0} such that IJagC ^(^) ^-^^ whose union over the elements 
doesn't cover all of L, we note that the innermost summation is run an even number of times 
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with the same term (there are 

2lAUaecMa)l equal terms). A gain, since the characteristic is 

two, these cancel. 

Qi9,r)= Ylg{a,h{a),r) 

Cecc{G) /i:C^2^\{0} clGC 

Restricted to the monomial for r'^', the only monomial (in r) in any g{a,X,r) which con- 
tributes to the final sum is the one for rl"^l, since there are no nonzero monomials of smaller 
degree in g{a,X,r), and the total degree of the monomial product should be L. We have 

[r^%i9, r)=Y E n Ka),r) 

Ea6clM«)l = |i| 

Since [r\^^]g{a, X,r) = f{a,X) and UaeC ^('^) ~ ^ together with ^aec\^('^)\ ~ l-^l i™plies 
\/a^b: h{a) n h{h) = 0, we get 

Cecc(G) h:C^2^\{0} a£C 
\/ai^b:h{a)r\h{b)=% 

which is the same expression as for [rl^l]p(/, r) at the end of the proof of Lemma [6l 

The above lemma offers an alternative route to compute [r^^^\p{f ,r) by evaluating the 
coefficient of r'^' in q{g,r) instead. To achieve this, we again compute q{g,r) for \L\n values 
on r and use Lagrange interpolation to retrieve the coefficient. To compute q{g, r) for a fixed 
r, we don't tabulate the values of g as we did for / in the evaluation of p{f,r). Instead 
we note that the inner sum in Y^Ycx'^g{-,-,r)0^)^ be evaluated in time and space 
polynomial in |A|. To see how, let Xi = A n V2 and X2 = A n and note that 

E = E + E ^9i.,.,r)iY) (7) 

YCX YCXi YCX2 

This decomposition is valid since g{a, Y, r) takes by definition the value for all Y such that 
both y n Ai 7^ and y n A2 / 0. The second sum in the rhs of ^ is easily evaluated since 
g{a, Y, r) also by definition takes the value zero for all Y such that \Y D A2I > 1. In the first 
sum of the rhs of ([7]) the coefficient of in row u and column v in the resulting matrix equals 
zero if u = V and otherwise evaluates to 

E E n ^'"^^ E n ^'^^ 

YCXi WeWn v{Y,l+l) wz&N W&^u wzi^W 

Sxiyv)=Y 

i.e. the number of weighted walks from n to u passing through / vertices (possibly with 
repetition) in X\ . These can be evaluated efficiently. Simple let an | X\ \ x | X\ \ matrix A with 
rows and columns representing vertices of Ai be defined by 

^ _{ Xuv ■■ uv e E 
^'^ \ : otherwise 
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and an |Vi| x \Xi\ matrix B with rows representing vertices of Vi and columns vertices of Xi 

by 

g ^ i Xuv ■■ uv G E 
'''^ 1^ : otherwise 

Then, the element at row u and column v of the matrix product BA'B^ by matrix multipli- 
cation definition equals the rhs of ([8]). This product can trivially be computed in time and 
space polynomial in n. 



5 TSP with Bounded Integer Weights 

In the TSP edges have weights I : E ^ and we seek the Hamiltonian cycle with the 
smallest total weight. Since our approach is algebraic, in particular operating on a sum- 
product ring, it is not evident how to handle a minimum query efficiently. For small weights 
though, there is an embedding solution. We note that the monomials in the polynomials we 
evaluate is a product over all edges along a Hamiltonian cycle. We introduce yet an auxiliary 
indeterminate y, meant to sort the Hamiltonian cycles after total weight. For every actual 
edge uv we represent, we also multiply it with y^('^^\ Formally, we extend the definition of / 
from Section 13.21 For uv £ F and C X C V2 , we set 

fy{uv,x)= Yl n y'^"'^^-^ 

For every arc uv £ F such that uv is an edge in G[Vi], and every d £ L^, we set 

/,(™,{d})=/(™)x,,,. 

Now X^™!^'"' A(Z), V2 U Lj, fy) seen as a polynomial in y, has a nonzero coefficient for the 
monomial only if there is a Hamiltonian cycle of total weight I. This idea of embedding 
a min-sum semi-ring on a sum-product ring is quite old. It was used by Kohn et al. [18] for 
TSP, and earlier by Yuval [23] for all-pairs shortest paths. 

We can use the Fast Fourier Transform on GF{2'^) to retrieve the smallest / for which 
the coefficient of is nonzero. First of all, we need to have 2^ larger than the maximum 
degree of a monomial in the indeterminate y to avoid having monomial coefficients wrap 
around We note that the largest degree equals the weight of the heaviest Hamiltonian 

cycle. This weight at least is less than w, the sum of all weights, and we set 2w > 2^ > w. 
For a generator g of the multiplicative group in GF{2^), we compute and tabulate T(/) = 
^m^ax U Li, fy) evaluated at y = for every / = 0, 1, • • • , 2'' - 2. The runtime of 

this step is in (9*(-u;2'^-^"'+™™°^). Next we compute the inverse Fourier transform of T: 

2* -2 

1=0 

The value t{j) equals the coefficient of y^ in XIS)'"" ^(^1 U Lj, fy), since for A; > 1 

1=0 ^ ^ 
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By the Fast Fourier Transform, t{j) is computed for all j in 0{wlogw) time. Finally, we 
search linearly in t{j) to find the smallest j such that the coefficient of is nonzero. 

To increase the probability of the lightest cycle to show in the sum (remembering that we 
only detect Hamiltonian cycles with at most rrimax arcs unlabeled by V2 this way), we rerun 
the algorithm r times with different partitions Vi U V2 = ^ as before. We output the smallest 
j found for which the coefficient of was nonzero in any run. This proves Theorem [3j 
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